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Abstract 

A generalized model taking into account the photoisomerization influence on 
the nematic ordering is presented. This generalized theory is used to investi- 
gate the effect of the nonextensivity on the concentration dependence of the 
long-range order parameter. The (/-dependent variation of the concentration 
of cis — trans isomers and of the order parameter with the time of exposure 
to the illumination is investigated within nonextensivity. It is also shown that 
for sufficiently long exposition of the mesophase to the illumination the ne- 
matic phase cannot disappear for some value of the entropic index, coming 
from Tsallis statistics. Moreover it is shown that long range interactions or 
the fractal structure in the liquid crystalline system might affect the charac- 
teristics of the phase transition in the physical system. We think therefore 
that this conclusion might shed light on the interaction potential energy terms 
in the similar system in future possible experiments, even to be performed on 
new objects with different symmetries. 
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I. INTRODUCTION 



A nonextensive entropy definition was proposed in 1988 [1] by Tsallis, 

1 _ V w n q 

S q = -k ^ Pl , qeR (1) 

where k is a positive constant, Pi is the probability of the system in the i th microstate, 
W is the total number of the configurations of the system and q is any real number. Since 
then, Tsallis thermostatistics (TT) has commonly been used for the investigation of the 
physical systems. As well-known, Boltzmann-Gibbs (BG) statistics is a powerful one to 
study a variety of the physical systems. However it falls for the systems which 1) have 
long-ranged interactions, 2) have long-ranged memory effects, and 3) evolve in multifractal 
space-time. In this manner, it has been understood that extensive BG statistics fails to study 
nonextensive physical systems not having these conditions. Consequently the standard BG 
thermostatistics is not universal and is appropriate for extensive systems. The parameter, 
q, in Eq.(l) is called the entropy index and measures the degree of the nonextensivity of 
the system under consideration. After a few years from 1988, the formalism was revised 
[2] introducing the unnormalized constraint to the internal energy. This generalization 
included the celebrated BG thermostatistics, which could be recovered when q — 1. After 
this work, this formalism has been commonly employed to study various physical systems. 
Some of them can be given such as; self gravitating systems [3,4], turbulence [5-8], anomalous 
diffusion [9-12], solar neutrinos [13], liquid crystals [14-18] etc. 

In spite of the successes of TT, some drawbacks were noted in the original formalism. 
These can be summarized as follows: i) the density operator was not invariant under a 
uniform translation of energy spectrum, ii) the g-expected value of the identity operator 
was not the unity, Hi) the energy of the physical system was not conserved. To overcome 
these difficulties, Tsallis et al. introduced the normalized constraint to the internal energy 
of the system [19]: 



The optimization of Tsallis entropy given by Eq.(l) according to this internal energy con- 
straint results in 
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where 



^ (l/Trm) (3 +(l-q)U^ P ^ 
Therefore the g-expectation value of any observable is defined as 



( Ai \ = b > ( 5 ) 

2^i=i Pi 



where A represents any observable quantity which commutes with Hamiltonian. As can be 
expected, the g-expectation value of the observable reduces to the conventional one when 
q = l. 

As for our motivation in this study, it was found in literature [20] that the structure of the 
mi cellar and liquid crystalline phases can be considered as self-similar. Based on this fact, 
recently a model of self-similar aggregation of alkylbenzenesulphonates (ABS) molecules 
was proposed [21]. Then the fractal structure of the sample was investigated experimentally 
[20] where the influence of the geometry of ABS molecules on the formation of their liquid 
crystalline phases was reported and the authors emphasized that ABS molecules are the 
fractal objects. As mentioned above, the nonextensivity may appear in systems where 
fractality exists. Therefore the main motivation of the present study comes from having 
fractal structure of some liquid crystals. 

As well-known from the literature and our last studies on nematic liquid crystals, Maier- 
Saupe theory (MST), which is a mean- field one, has some insufficiencies coming from; 1) as 
being a mean-field theory, 2) neglecting higher interaction terms in the interaction potential 



function. Due to these insufficiencies, MST gives, for instance, the critical value of the 
long range order parameter P 2 = 0.429 for all nematics. Considering these terms (i.e. 
the possibility of the fractality in liquid crystals and the insufficiencies of the mean-field 
approach), we apply TT to the nematic liquid crystals. To this end, in [14], MST was 
generalized within TT and it was showed the influences of the nonextensivity on the nematic- 
isotropic transition. Next the dimerization process was handled using TT in [15] which can 
be regarded as the first application of TT to a liquid crystalline system. In [15], the influence 
of the dimerization process on the nematic ordering is investigated by using a nonextensive 
thermostatistics (TT). A theoretical model taking into account the dimerization influence on 
the nematic scalar order parameter has been summarized and theoretical predictions for the 
nematic order parameter are improved by using TT. Then the P 4 model with P 4 -interaction 
(P 4 indicates the Legendre polynomial for L — 4) was investigated and generalized within 
TT in [16]. The dependence of the second order parameter (P 2 ) on the fourth one (P 4 ) 
was studied in [17] and the generalized theory was applied to PAA (para-azoxyanisole), a 
nematic liquid crystal. In [18], a generalized mean-field theory was presented relating helix 
tilt in a bilayer to lipid disorder. It has been observed from all these studies that if one uses 
a generalized form of the MST within TT, the deviations between the theory (MST) and 
the experimental data are minimized and the obtained results agree with the experimantal 
data very well. 

In this study, the influence of the photoisomerization on the orientational ordering is 
handled by using a generalized theory. The standard theory has been recently proposed in 
literature [22] . It is an extension of the MST and could be an appropriate one to investigate 
the variation of the concentration of cis — trans isomers and of the scalar order parameter 
with the time of exposure to the illumination. As expected, as it is an extension of a 
mean-field theory (MST), it might have some deviations form the experimental data; for 
instance it also gives a universal value of the long range order parameter at transition point, 
0.429. To overcome similar insufficiencies (considering the possible fractality of the liquid 
crystalline system), we use a generalized form of the standard theory and investigate the 
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physical quantities concerning physical system. To this end, we first summarize the standard 
and generalized theories. Then we investigate the effects of nonextensivity on the nematic- 
isotropic transition and also on the evolution of the concentration of cis — trans isomers and 
of the scalar order parameter with the time of exposure to the illumination. 



The standard theory [22] assumes that the liquid crystalline system is formed by a 
mixture of two elements, namely two miscible nematogens Aq and N 2 . According to the 
system under consideration can be dimers and monomers, if an energy to break the hydrogen 
bonding is considered [23], or trans and cis isomers in the case of azobenzene materials [24]. 
Both kinds of molecules, namely Aq and N 2 , are supposed to be rodlike. In the mean-field 
approximation, the total nematic potential function can be given by [25-27] 



where 9\ and 9 2 are the angles formed by the molecular long axes of the molecules of Aq 
and A*2 with the nematic director of the mixture respectively, Vu and Vi 2 are the mean- 
field potentials acting on a molcule of Aq due to the other molecules of Aq and due to the 
molecules of N 2 respectively. A similar meaning has V 21 and V 22 . These potential functions 
can be written as 



in the MST, where a^'s are coupling constants and they depend on the distance between 
the center of mass of the molecules, on a molecular property and also on the concentration 
of one of the constituents in the system [23]. Using these terms above, we can write the 
nematic potential function as 



II. THE GENERALIZED THEORY 



V (0i, 2 ) = XX- = v u + v l2 + v 21 + v 2 



22, 



(6) 



Vy = -a t3 P 2 (9,) (P 2 (9,)) , 



(7) 



V(9 1 ,9 2 ) =V 1 (9 1 ) + V 2 (9 2 ) 




(8) 
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where P 2 (9i) = (^(^i)) i s the scalar order parameters of the molecules of JVj. The scalar 
order parameters are determined by solving the set of coupled equations according to Eq. (5) 
for i — 1,2 self consistently. These coupled equations depend on the temperature T and on 
the quantities appearing in the nematic potential function. The quantities contain 
the strengths of mean- field interaction among the molecules forming the nematic 1, among 
the molecules forming the nematic 2, and among the molecules of 1 and those of 2. a'^s are 
assumed as 

OCij UijTlj^ (9) 

where tty = Uy t and rii = N i /{N 1 + N 2 ), with iVj being the number of molecules of the 
nematic % [23]. 

The generalized free energy of the system is given by 

F i = ~ E (\ N i M + N * kT M^)) , (10) 

where ln q (Z iq ) = (Z}~ q — 1)/(1 — q) and (T^) = —J2j a ijP2(Qj)P2(6i)- Since we assume 
uu = u 2 i = u, the free energy density of the system can be given by 



F * = N { - E [l n l u a ( p 2 (0 t ) f + n t kT h^iZi) 



+ un l n 2 P 2 {6 1 )P 2 {6 2 )\ , (11) 



where N = N x + N 2 . 

Let us now consider the case in which one of the constituents do not contribute to the 
nematic order. The experimental evidence shows that the molecules of the cis isomers are of 
banana shape. According to recent investigations, this banana shape can be approximated 
by two linear parts forming an angle between them not very far from n/2 [22]. As a result 
of this configuration, the molecules of the cis isomer behave as spherical objects due to the 
spinning and tumbling motions of thermal origin. Consequently in the liquid phase they 
can originate just an isotropic phase. In this manner, the molecules of cis isomer lead to a 
kind of neutral background, where the molecules of trans isomer are responsible for nematic 
order. So we will consider that the trans isomers have cylindrical symmetry and can exhibit 
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a nematic order, whereas the cis isomers will tend to reduce this orientational order. In 
addition, we also assume that at zero temperature the system under consideration is made 
by Ntrans = N isomers. At a given temperature T we have 

^trans ~l" n cis lj (12) 

where n tra ns and n cis are the equilibrium concentration of trans and cis isomers and are 
given by 

^trans trans \I) /-A 

n cis = N as (T)/N 

at a given temperature respectively. Thus the quantities are reduced to the case of trans 
isomer-trans isomer interaction, i.e., = a tra ns and a trans = un trans [23]. 

In the case of one component, P2(#i) = (-fM#)) will be the scalar order parameter of 
trans isomers. The standard theory is based on MST with a ingredient represented by 
temperature-dependent concentration. This dependence in the standard theory plays an 
important role and it departs from the one represented by the concentration dependence 
characteristic of the excluded volume theories of Onsager type [28,29]. As well-known, in 
Onsager model the number of particle is fixed. In this manner, the critical density, below 
which the isotropic phase is stable, follows from a balance between the energy connected 
with the excluded volume and the thermal energy. In contrast, according to the standard 
model which will be generalized below, the number of particles responsible for the nematic 
phase depends on the temperature. In Onsager's theory, this is equivalent to having an 
excluded volume that is temperature dependent. 

In Fig.(l), we show the dependence of the long range order parameter on the concen- 
tration of cis isomers for a fixed reduced temperature T R = 0.1 u/k and for some q values. 
The nematic scalar order parameter vs. the concentration of the cis isomers in the system 
for Tr = O.lu/k (reduced temperature) is illustrated in Fig.(l). Since the standard theory 
(q = 1) depends on Maier-Saupe mean-field interaction energy, the system, as expected, ex- 
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hibits a first-order phase transition to the isotropic phase for n c cis = 0.55 and \^P 2 j decreases 
suddenly to zero from 0.429 which MST gives for all nematic liquid crystals. However it 
is experimentally known [30] that the critical value of the nematic order parameter varies 
from one nematic to the other and it lies in a range of 0.25 — 0.5 for nematic liquid crystals. 
We think therefore that the standard model based on mean-field approach appears to be 
insufficient and alternative approach is needed. On the other hand, as q is varied, we observe 
that the critical value of the order parameter changes. The obtained results are summarized 
in Table in which the entropy change at transition point is also reported as a function of q. 
We denote the entropy of the system as S q . 

In Fig. (2), the phase diagram of temperature vs. concentration of cis component for the 
system is plotted. One observes that the usual critical temperature (T R ~ 0.2202-u/A;) for 
MST [31] is reached when n cis = as expected, since the concentration of the molecules 
leading to the nematic phase is fixed in MST. An interesting point is that the experimental 
data reported in [25] show that the critical temperature for the nematic-isotropic phase 
transition of the mixture is a linear function of the concentration of the solute. In this 
framework, we observe from Fig. (2) that the entropic index q just affects the slope of the 
line associated with this linear dependency. 



III. ORDER PARAMETER VS. EXPOSURE TIME FOR ILLUMINATION 

At a given temperature, the isomerization reaction follows the reaction scheme trans — > 
cis, depending on some rate constants that, in turn, depend on the intensity of the light in the 
illumination process [32,33]. Then time variation of the concentrations of two constituents 
can be approximated by 

(T,t)=cn trans (T) + (l-c) 

^trans 

(T) exp^H/r) (13) 

within TT, where r is a characteristic time and c is a parameter to control the fraction 
of trans-cis isomers at a given temperature after the illumination. When t — 0, we have 
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n tra ns(T, 0) = n tran s{T). From ntrans + n cis = 1 and assuming for n tra ns the approximate 
expression given by Eq.(13), we obtain 

U cis (T, t) = 1 - Utrans (T, t). (14) 

It is the experimental fact that with the exposure time to the illumination, the concentration 
of cis isomers increases, whereas the trans isomer concentration does not tend to zero, but 
to some fraction, controlled by the parameter c, of the equilibrium concentration n trans (T). 

Now we would like to discuss about the equilibrium distribution of the concentrations 
of the different constituents forming the mixture. At a given temperature, the equilibrium 
distributions of cis and trans isomers in TT could be given by 

n cis = exp q (-P(fi + E)) 
n trans = exp q (-P(i) (15) 

respectively, where \i is the chemical potential of the mixture. Now if the above results 
are used in Eq.(12), it is possible to observe the behaviour of the order parameter with 
the exposure time to the illumination. In Fig. (3), we illustrate the variation of the order 
parameter with the exposure time for the reduced temperature of Tr = O.lu/k, c = 0.3 and 
for some q values. We observe that for a given value of c, the entropic index q plays an 
important role and has a considerable effect on the variation of scalar order parameter with 
the exposure time. The similar behaviour can be easily seen in Fig. (4) in which the order 
parameter vs. exposure time is plotted for c = 0.4 and some values of q. Probably the most 
interesting result is seen in Fig. (5); for q — 1 (standard theory) and c > 0.55, there is no 
phase transition at temperature T R = O.lu/k. However for q < 0.5 and c = 0.55, a first 
order phase transition occurs at this temperature (Tr = O.lu/k). This conclusion implies 
that the permenance of the liquid crystalline phase depends not only on the concentration 
of cis isomers in the system, but also on the entropic index q (i.e. probably long range 
interactions or fractality). 
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Q 


0.98 


0.99 


1 


1.01 


1.02 


Sc 


0.346 


0.391 


0.429 


0.463 


0.49 


AS q /Nk 


-0.281 


-0.353 


-0.418 


-0.479 


-0.528 



Table. 
IV. SUMMING UP 



We have presented a generalized model to study on the behaviour of the scalar order 
parameter in liquid crystalline systems, consisting of more than one nematogen component. 
It is the experimental fact that some liquid crystals have fractal structure. This fact implies 
that TT can be applied to the liquid crystalline systems. In addition, the standard model, 
which has been generalized, is an extension of MST and thus some drawbacks coming from 
insufficiencies of mean-field approach could be expected. In this framework, a generalized 
model which is an extension of MST has been used. The generalized theory can be success- 
fully applied to distinct systems, such as one formed by cyclic dimers and monomers, and 
also to the systems where azobenzene compounds play an important role. For example, a 
similar system, consisting of closed and open dimers, has been studied in [15], where the 
experimental data related to the dimeric system has been successfully explained by using a 
generalized theory. It must be noted that the success of the generalized mean-field theories 
in our earlier studies may come from being fractal structure of liquid crystals, apart from 
considering long range interactions. Similar successful explanations can be expected in other 
dimeric and liquid crystalline systems. 

It has been also studied on the systems where the illumination process gives rise to the 
appearance of new objects with different symmetries in the systems. We have shown that 
fractality or long range interactions (i.e. q ^ 1) might affect the characteristics of the phase 
transition in the physical system. This conclusion might shed light on the interaction poten- 
tial energy terms or fractal structure in the similar system in future possible experiments, 
even to be performed on new objects with different symmetries. 
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FIGURE AND TABLE CAPTIONS 



Figure 1. The variation of the scalar order parameter with the concentration of cis isomer 
for a fixed reduced temperature T R = O.lu/k and for some q values. As q increases, the 
critical value of the scalar order parameter at transition point also increases. 

Figure 2. g-dependent phase diagram in the reduced temperature vs. cis concentra- 
tion plane. The critical lines represent the first order phase transition, namely nematic- 
isotropic transition, induced by the variation in the concentration of cis isomers in the 
nematic medium. 

Figure 3. The variation of the scalar order parameter with the exposure time for uv 
illumination for c = 0.3 and some q values, r is a characteristic time. As seen from figure, 
the entropic index q affects the variation of scalar order parameter with the exposure time 
considerably. 

Figure 4. The variation of the scalar order parameter with the exposure time for uv 
illumination for c = 0.4 and some q values. As seen from figure, the entropic index q affects 
the variation of scalar order parameter with the exposure time considerably. 

Figure 5. The variation of the scalar order parameter with the exposure time for uv 
illumination for c = 0.55 and some q values. As seen from figure, whereas there is no phase 
transition in the standard theory (g = 1) for c > 0.55, the generalized theory assumes a first 
order phase transition for q < 0.5. 

Table. Entropy change of the system at transition point for some values of the entropic 
index q. 
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